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Внедрение кругового штампа при заданной постоянной температуре 
на плоской подошве штампа в непрерывно неоднородное 
термоупругое полупространство” 


Л. И. Кренёв, С. М. Айзикович, Б. И. Митрин 


Рассматривается осесимметричная квазистатическая задача термоупругости о внедрении цилиндрического 
штампа с плоской подошвой, на которой поддерживается постоянная температура, в функционально- 
градиентное полупространство, модуль упругости, коэффициент Пуассона, коэффициенты теплопроводности 
и линейного расширения которого непрерывно изменяются в приповерхностном слое независимо друг от 
друга. Вне контактной зоны поверхность идеально теплоизолирована и свободна от напряжений. При реше- 
нии задачи используются полученное ранее с помощью численно-аналитических методов (аппарата инте- 
гральных преобразований Ханкеля и метода модулирующих функций) решение несмешанной задачи о про- 
извольном термомеханическом воздействии на неоднородное по глубине термоупругое полупространство. 
Решение исходной задачи сводится к решению системы парных интегральных уравнений. Свойства транс- 
формант ядер парных интегральных уравнений задачи позволяют применить хорошо обоснованный и разви- 
ваемый в настоящее время двусторонний асимптотический метод. С помощью данного метода найдены в 
аналитическом виде приближённые выражения для величин теплового потока и смещения поверхности по- 
лупространства, контактные напряжения под подошвой разогретого штампа. Приведены значения контакт- 
ных напряжений для различных случаев изменения механических и температурных свойств в приповерх- 
ностном слое. Рассматриваются случаи, когда значения термоупругих свойств покрытия совпадают со значе- 
ниями термоупругих свойств подложки, либо когда значение характеристики отличается в 2 раза (в большую 
или в меньшую сторону) на поверхности и линейно убывает (или растёт) по глубине до значения в подложке. 
Ключевые слова: смешанные задачи, неоднородные материалы, термоупругость, функционально-гради- 
ентные материалы, аналитические методы. 


Введение. Учёт неоднородности свойств материала при моделировании тепломеханического 
воздействия является актуальной задачей термоупругости. Такая задача возникает, например, 
при исследовании свойств защитных покрытий из функционально-градиентных материалов для 
частей машин и механизмов, подверженных интенсивному термомеханическому воздействию. 

Первые аналитические решения задачи об индентировании полубесконечного изотропно- 
го упругого тела осесимметричным нагреваемым штампом были получены в работах [1, 2]. Свя- 
занным задачам термоупругости были посвящены монографии В. Новацкого [3] и В. Г. Карнаухова 
[4]. Контактные задачи термоупругости рассматривались в статье [5] и последующих работах 
Д. В. Грилицкого и его соавторов. В работах Дж. Барбера (например, [6]) изучались задачи, свя- 
занные с контактом движущихся частей механизмов. Свой вклад в развитие термоупругих кон- 
тактных задач также внесли К. А. Вицоп, М. Мода, В. Паук и многие другие. Из последних работ, 
рассматривающих задачу о внедрении горячего штампа, можно отметить [7]. В этой работе полу- 
чено аналитическое решение задачи для трансверсально изотропного полупространства и рас- 
смотрено его применение для сканирующей зондовой микроскопии. 

В перечисленных выше работах предполагалось, что свойства деформируемой среды од- 
нородны по её объёму, и, таким образом, они не всегда являются применимыми для случая де- 
формирования тел с покрытиями или функционально-градиентных тел. При этом современные 
экспериментальные исследования показывают эффективность их применения в промышлености. 


* Работа выполнена при финансовой поддержке Министерства образования и науки РФ (соглашение № 14.В37.21.1632) и 
РФФИ (гранты № 12-07-00639-а, № 13-08-01435-а). 
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Например, 2. О. И в работе [8] рассмотрел температурные свойства \57-№Сг с промежуточным 
слоем, состоящим из ФГМ; им была показана высокая температурная устойчивость материала. 

Интерес к моделированию термоупругого поведения ФГМ существовал с момента начала 
их использования в промышленности. В последнее десятилетие учёными по всему миру был по- 
лучен ряд результатов [9-10], в которых различными методами строились математические моде- 
ли такого поведения при различного вида нагрузках. Однако решения, используемые в приведён- 
ных работах, используют различные упрощённые модели изменения характеристик в объёме ма- 
териала, и не могут быть применены для сложных случаев изменения термоупругих и механиче- 
ских свойств по глубине слоя. 

В настоящей работе рассматривается осесимметричная квазистатическая задача термо- 

упругости для функционально-градиентного полупространства, модуль упругости и коэффициен- 
ты Пуассона, теплопроводности и линейного расширения которого непрерывно изменяются в 
приповерхностном слое независимо друг от друга. Для её решения используется численно- 
аналитический метод, описанный в работе [11]. При этом используется полученное ранее [12] 
решение для определения искривления формы поверхности неоднородного по глубине термо- 
упругого полупространства при локальном нагреве. 
1. Постановка задачи о воздействии горячего штампа на термоупругое неоднородное 
покрытие, сцепленное с полупространством. Рассмотрим полупространство ©, термомеха- 
нические характеристики которого непрерывно меняются с глубиной в пределах прилегающего к 
поверхности слоя толщины Н, а затем стабилизируются и остаются постоянными. С полупро- 
странством свяжем цилиндрическую систему координат (г, ф,27). При этом полагаем, что 
О<г<о, а 2<0. Обозначим через и‚,и,и/ смещения вдоль осей г,ф,2, а через 
0,0, 0б,,Т.,Т,, То, — Радиальное, угловое, нормальное и тангенциальные напряжения соот- 
ветственно. Кроме того, введём обозначения для температуры — Т, коэффициента теплопро- 
водности — Л, (2), коэффициента линейного расширения — а; (2), коэффициента теплоёмко- 
СТИ —*:6- (2). В ненапряжённом состоянии температуру полупространства примем равной Т, и 
будем анализировать разность температур 9 =Т -Т.. 


На поверхность Г неоднородного полупространства © воздействует жёсткий круговой в 
плане штамп-теплоизолятор, который вдавливается в полупространство силой Р. На контактной 
площадке 0<г<а поддерживается постоянная температура и образуется тепловой поток, 
направленный внутрь неоднородного полупространства. 

Далее будем полагать, что коэффициенты Ламе М(2) и Л(2), и коэффициенты тепло- 


проводности Л, (2) и линейного расширения а, (2) являются непрерывными функциями коорди- 


наты 2 , такими, что 


2.М(2) =М (2), Л(2)=^° (2), -Н<2<0 
аЕ Ааа, (28а: (2). 
3.М° (-Н) = М, ^° (-Н) = Л, 
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т Л(2)>А, >0, пах ^(2 


2Е(-о, 0) 2«( —, 0) )= Л 
пи М(2)>М, >0, пах М(2)<М 
я 


2Е(-о, 0) 2«( —, 0) 


(1) 
пи Л, (2)>Л,. >0, „тах, Л; (2)<^; 


2Е(-о, 0) 2Е(-о, 0) 


„ПМП б; (=) а. =0; „пах, а; Е а 


где Н — толщина неоднородного слоя, сцепленного с подстилающим полупространством, то есть 

глубина, с которой мы полагаем термомеханические характеристики полупространства посто- 

янными; индекс $ соответствует подстилающему однородному полупространству; С — не- 
однородному слою; Л., Л", М,, М", Л;., Л;,а;., а; — произвольные константы. 

Наряду с парой коэффициентов Ламе для описания упругого поведения твёрдого изотроп- 

ного тела используются: модуль сдвига с и коэффициент Пуассона у или модуль Юнга Е и ко- 

эффициент Пуассона у. Коэффициенты Ламе Л и М (иногда обозначаемый С и называемый мо- 


дулем сдвига) связаны с модулем Юнга Е, коэффициентом Пуассона \ и изотермическим моду- 
лем объёмного расширения К соотношениями 


_ Е й ЕУ _ М(2М+3Л) 
_ 2(1+%)° (1+\)(1-2%)' М+Лл ю 
“= мел), К (2) =3^(2)+2м(2). 


Под действием центрально приложенной силы Р штамп переместится в направлении оси 2 на 
величину д. Граничные условия при сделанных предположениях имеют вид: 


0. =0, гла, 
и=-5(г)=-(5-4(г)), г<а, 

2=0, т, =Ти =0, 10 =0,, г<а, (3) 
И гра. 
02 


На границе сцепления неоднородного слоя с однородным полупространством, при 2 =-Н в силу 
непрерывности, должны выполняться условия сопряжения по смещениям, их производным, тем- 
пературе и тепловому потоку. 
С 5 С 5 С 5 
0° (г.-Н)=с (г,.=Н), т (.-Н)=т» (г.-Н), и (г,-Н) =и (г.-Н), 


2 2 


ше (,-Н)=и? (г,-Н), 60° (г,-Н)=65 (^,-Н), (6) (^,-Н) = (6%) (=,-Н). 


На бесконечности при (г, -2)>® значения величин смещений, деформаций и напряжений 


(4) 


стремятся к нулю. Значения разности температур и теплового потока при этом также стремятся к 
нулю. 
ПР (ОИ бб быта 0,0) =0, 


А. (5) 
Ни (М.И, =, еб О т 9, 6') =0. 


2—0 
Требуется определить распределение контактных нормальных напряжений под штампом и вели- 
чину теплового потока в пределах контактной зоны: 


171 777 





#0 =-Р(г), аа. 
6 
22|. (0), г<а (6) 
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Рис. 1. Схема внедрения горячего штампа с плоской подошвой в неоднородное полупространство 


2. Сведение задачи о внедрении гладкого кругового горячего штампа с плоской по- 
дошвой в неоднородное полупространство, к решению системы парных интегральных 
уравнений. При построении системы парных интегральных уравнений осесимметричной сме- 
шанной задачи статической термоупругости воспользуемся полученными в работе [12] формула- 
ми, выражающими зависимость величины смещения и/ и температуры 0 через неизвестное рас- 
пределение контактного давления р(г) и теплового потока через контактную зону 
00(г, 2) 
и 


2=0 


и (0) => 15 м урар[ (у, 0) 2 (\р)-» (чг)ау 


[=1 


+9 в (Ур) 2. (гам + [В(о)рав[ ву, 0) 2, (Ур) 2, (уг) ау, (7) 


6(^,0)= 15 (г,0)= [В() ии 


0 
Так как мы полагаем, что контакт штампа с неоднородным основанием гладкий, то отсюда следу- 
ет равенство нулю значений 9(г) =0, О<г<а, атакже второго интеграла в выражении (7). Те- 
перь условие (7) можно записать в виде: 


и (г, > )= (р зрар[ (У, 0) 2 (Ур) 7% (уг) ау + 


+[В ты | (у, 0) 2, (ур) 7. (уг)ау, (8) 


0 


9(-,0) = 1. (г, 0) = 


5—ъ-\ 


В(р)рар[ в (У, 0) 1 (уг) ау. 


В рассматриваемой нами задаче упругие характеристики не зависят от температуры, поэтому ре- 
шение системы (8) разбивается на два этапа. 

На первом (рассмотренном в работе [12]) определяется величина теплового потока и рас- 
пределение температуры в неоднородном полупространстве, а также деформация поверхности 
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неоднородного термоупругого полупространства при его нагреве с поверхности источником с по- 
стоянной температурой. 

В настоящей работе решается смешанная задача о внедрении гладкого кругового штампа 
с плоской подошвой в неоднородное полупространство, деформированное воздействием найден- 
ного нами теплового потока через контактную зону. 

Заметим также, что в связи с отсутствием сцепления индентора с поверхностью полупро- 
странства, для корректной постановки задачи должно быть выполнено соотношение р(а)>0 для 


всех г<а. 


Для решения поставленной задачи приведём её к эквивалентному парному интегральному 
уравнению относительно р(г), для этого воспользуемся результатами работы [11]. Запишем 


граничное условие на поверхности слоя в пределах контактной зоны, при этом обозначим 
[, (и, 0) через [ (и): 





ИА 0551, 
(9) 
[Р(о) (аг) ада =0, ПЫ[, 
}- [6 ‚ (ар)рар, НО 


и 5=х/а, 0<г<1. 
Здесь Л, — функция Бесселя нулевого порядка, ф(га) =0 для плоского штампа, 9 — величина 
внедрения штампа в полупространство, Л =Н / а — характерный геометрический параметр зада- 
чи, [ (Ла) — трансформанта ядра парного интегрального уравнения. При выполнении условий: 


тм Е (2)(1 ты >С, >0, 


2е(-о; 0] 


тах Е (2)(1-\)"<с, <®, 
2«(-;0] (2 )( ) 2 (11) 
[0 Е(2)(1-У)` = СОП$. 


р——0 


Аналогично [15] может быть показано, что трансформанта ядра [ (а) обладает следующими 


свойствами: 
[(а)= А-В, а! +О(а*), а 0, 
К-во) ы 
Е (а) =1+В, а] +0(а?), а-> <, 
А= Шт Е (0)Е (2) `. (13) 


Здесь В,,В, — постоянные. Свойства (12)-(13) функции [(а) позволяют аппроксимировать её 


дробно-рациональным выражением специального вида. 
Введём следующие определения. Будем говорить, что (У) принадлежит классу П)„, , если 


она имеет вид: 
№ = 
в, (у) = (№ +^)(\ +в?) ”, (в,-В,)(А-А,)=0, =. (14) 
1=1 
Здесь А,В, (1=1,2,..., М) — некоторые комплексные постоянные. 
В свою очередь будем считать функцию [ (У) принадлежащей классу >„, если её можно 


представить в виде: 
М 


1: (у) => Си (у +22]. (15) 


К=1 
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Здесь С, — действительные константы, а р, (К =1,2,..., М) — комплексные. 

В работах [13, 14] доказано, что если функция [ (а) обладает свойствами (12)-—(13), то 
она может быть представлена в виде: 

1 (а^) = [и (а^) + 1» (а^). (16) 

3. Построение приближённого аналитического решения задачи о внедрении нагретого 
штампа в термоупругое полупространство. Решение (9) ищется следующим образом: парное 
интегральное уравнение приводится к более простому виду, ядро аппроксимируется дробно- 
рациональным выражением, с помощью операционного исчисления находится решение, а затем 
обратными преобразованиями задача сводится к решению линейной алгебраической системы. 

Рассмотрим парное интегральное уравнение вида: 


ГР (У) 2, (гу) ау = 6,5(1+! ()),г <1, 
(17) 
ГР(У) о (ГУ) )уау = 0, РЕ 


Функция Г (г) описывает деформацию поверхности под действием источника тепла с постоянной 
температурой. Её можно представить в виде ряда Фурье-Бесселя: 


г (г) => 6.2» (г). (18) 


Учитывая, что при = > 0, 7, (/=) > 1, уравнение (17) посредством операторов: 


Уф(г еее г; (9:2, (а = с0$ аб), 
уф (г) ча Е-т®а г; (95. (аг) = а* созаб), 


можно преобразовать к виду: 


а5[Р(9 1 (Лу) созёуау = У, со, Ё, 0 <Ё<1, 





(19) 
[Р(у)соз уду =0, Ё>1. 
0 
В случае, когда 
В (№) м му? + А? 
Е(Лу) =Ё (Л ; 20 
(ЛУ) » (Ау) = Р; (у *) И вут+ в ( ) 
можно получить решение задачи в аналитическом виде. 
Введём функцию 
р(#) = [Р (а) созаа. (21) 
0 


Тогда (19) можно представить в виде: 


а? 
(5) р(#)= 2, (2), (#),р= Ч, Е =[0;1], (22) 
где полиномы Р и Р, определяются, как в выражении (20). 
Решая неоднородное дифференциальное уравнение (22) и выполняя затем обратное преоб- 
разование Фурье для (21) получим выражение для контактных напряжений под горячим штампом: 
М со 
29 УСА А+ ЬЖ (ии Е (м (23) 
п 7 71 


г 
39 





В 1-г 
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Здесь введены обозначения: 
_ А __ | Е с0$| т эп СЕ 
Ф (г, А) = ГМ) = | 
(А) Ты Ая Е" о 1 - г? ео 
Постоянные С, определяются из системы линейных алгебраических уравнений: 


са, ыы ь В ‚5 0, К=1,2,...,М, (24) 





где 
ВСПА+ АВА. В со5и- ияйр 
ОВ а РВ а 
В“-А [ (Ан) (ны + В ) 
В данном случае связь между приложенной силой и осадкой штампа имеет вид: 
р-4е, Г (0 ) + СА“ ПА" +В) Г (^н ‚)ну пн, у (25) 


4. Определение контактных напряжений под поверхностью горячего штампа, внедря- 
емого в неоднородное термоупругое покрытие. Используя выражение (23), рассчитаем зна- 
чения контактных давлений под горячим штампом для 5 характерных случаев изменения свойств 
по глубине: 
0 — значение характеристики в покрытии не отличаются от соответствующих значений в подложке; 
3 — значение характеристики на поверхности в 2 раза превышает значение в подложке и линей- 
но убывает по всей толщине покрытия до значения в подложке; 
4 — значение характеристики на поверхности в 2 раза меньше значения в подложке и линейно 
возрастает по всей толщине покрытия до значения в подложке. 

Наглядное изображение законов изменения термомеханических характеристик приводится 
на рисунке 2. 
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Рис. 2. Схема характерных случаев изменения свойств в термоупругом покрытии 


На рис. 3 представлены значения контактных напряжений под поверхностью горячего 
штампа с плоской подошвой, внедряемого в непрерывно неоднородное полупространство. Обо- 
значения кривых соответствуют законам изменения термоупругих характеристик, описанным вы- 
ше. Максимальный уровень напряжений достигается в том случае, когда модуль Юнга и коэффи- 
циент Пуассона возрастают по глубине покрытия, а коэффициент теплопроводности и коэффици- 
ент линейного расширения убывают. 
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Рис. 3. Контактные напряжения под поверхностью горячего штампа, внедряемого в неоднородное 
термоупругое полупространство 


С учётом независимого изменения модуля Юнга, коэффициента Пуассона, коэффициента 
теплопроводности и коэффициента линейного расширения в покрытии, при рассматриваемых 
нами законах линейного убывания и возрастания в 2 раза относительно свойств подложки число 
различных комбинаций составляет 16. Кроме этого, рассмотрим случай однородного полупро- 
странства. На рис. 4, 5 приведены отношения контактных напряжений на поверхности к контакт- 
ным напряжениям для однородного полупространства в случае внедрения в него горячего штам- 
па с плоской подошвой для 3 отношений радиуса зоны контакта к толщине неоднородного слоя. 
Обозначения экземпляров наборов независимого изменения неоднородных термоупругих харак- 
теристик на рисунках описаны выше. 
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Рис. 4. Относительные контактные напряжения под поверхностью горячего штампа, внедряемого в неоднородное 
термоупругое полупространство. Отношение радиуса зоны контакта к толщине слоя равно 0,25 
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Рис. 5. Относительные контактные напряжения под поверхностью горячего штампа, внедряемого в неоднородное 
термоупругое полупространство. Отношение радиуса зоны контакта к толщине слоя равно 4,0 


Анализируя рисунки можно отметить, что в случае покрытия, более жёсткого на поверх- 
ности, контактные напряжения в центре зоны контакта меньше, чем в однородном полупростран- 
стве, а для мягких на поверхности покрытий — больше. С увеличением зоны контакта разброс 
кривых увеличивается. 

Заключение. С помощью двусторонне-асимптотического метода получено эффективное анали- 
тическое решение осесимметричной контактной задачи о внедрении кругового штампа с плоской 
подошвой при заданной на ней постоянной температуре в непрерывно неоднородное термоупру- 
гое полупространство. Были получены численные результаты для относительных контактных 
напряжений, которые позволяют судить о влиянии характера изменения упругих модулей и тер- 
момеханических коэффициентов в покрытии на величину контактных напряжений. Данные могут 
быть использованы при разработке функционально-градиентных покрытий материалов, исполь- 
зующихся в условиях интенсивного термомеханического нагружения. 
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ЕЕАТ ВОТТОМ^ 


Е. Г. Кгепет, $. М. АхШКомсН, В. Т. Митт 
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{ептотесвапга! ргорегйез уапаноп т {те Боипаагу !вуег оЁ #1е Ва!5расе аге ргоиЕд. ТПе сазез ейПег иеп 
уа/иез о! {ве Пептотесвагтгса! соа па ргорег@ез аге {пе 5ате аз {Позже оГЁ те 5ирзёгее, ог ийеп {те ргорейу 
уа/ие аТег$ виисе (ириага ог аоигпигага) оп пе 5итасе, апа пеапу аесгеазез (ог доез ир) т ЧерЁ ю те уаие 
т ше из, аге соп'аегед. 

Кеуигога5: ти/хеа ргоМетз, пвотодепеоиз тсйепаБ, Пегтое/азйсйу, ГипсИопа/у-дгадея тейепаБ, апг!уйс 
те#осе. 


* Тпе гезеагсп {5 допе м ИИ {Не Япапаа! зирро! #гот Не Ми гу оЁ ЕдисаНоп апа З&епсе оЁ {пе Киз$ап ЕедегаЧоп (Адгеетеги 
по. 14.В37.21.1632), апа от ВЕЕТ (Сгап$ поз. 12-07-00639-а, 13-08-01435-а). 
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